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Orchestration of signaling, photoreceptor structural integrity, and maintenance needed for 
mammalian vision remain enigmatic. By integrating three proteomic data sets, literature mining, 
computational analyses, and structural information, we have generated a multiscale signal 
transduction network linked to the visual G protein-coupled receptor (GPCR) rhodopsin, the major 
protein component of rod outer segments. This network was complemented by domain 
decomposition of protein-protein interactions and then qualified for mutually exclusive or 
mutually compatible interactions and ternary complex formation using structural data. The 
resulting information not only offers a comprehensive view of signal transduction induced by this 
GPCR but also suggests novel signaling routes to cytoskeleton dynamics and vesicular trafficking, 
predicting an important level of regulation through small GTPases. Further, it demonstrates a 
specific disease susceptibility of the core visual pathway due to the uniqueness of its components 
present mainly in the eye. As a comprehensive multiscale network, it can serve as a basis to 
elucidate the physiological principles of photoreceptor function, identify potential disease- 
associated genes and proteins, and guide the development of therapies that target specific branches 
of the signaling pathway. 
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Introduction 

The work of many different groups over the past decades has 
led to a detailed understanding of the molecular mechanisms 
underlying the initial steps of the vision process in photo- 
receptor cells (Palczewski, 2006; Kwok et al, 2008; reviewed in 
Ridge et al, 2003) . Rod photoreceptor cells are neurons capable 
of converting light into electrical signals. They possess a 
specialized structure consisting of five principal regions 
(Figure 1A): (i) the rod outer segment (ROS) composed of 
~800 closed membrane discs where phototransduction takes 
place; (ii) the connecting cilium (CC) that joins the outer 
segment to the rest of the cell and regulates the traffic of 
proteins and other components in both directions; (iii) the rod 
inner segment (RIS) responsible for general cell metabolism, 
housekeeping, and protein production; (iv) the cell body with 
the nucleus (N) ; and (v) the synaptic region (SR) that makes 
the electrical connections to the neurons in the retina. Protein 



activity and turnover in the ROS are highly dynamic: about 
10% of all discs are generated each day at the base of the 
segment, while older discs are removed at the distal end by 
phagocytosis of the neighboring retinal pigment epithelium 
cells (Boesze-Battaglia and Goldberg, 2002) . To replenish, the 
components of the ROS and the vesicles synthesized in the RIS 
compartment need to be transported through the CC region, 
either actively or by diffusion (Reidel et al, 2008). 

Rhodopsin is the major visual pigment in rod photoreceptor 
cells. It is a prototypical seven transmembrane-spanning G 
protein-coupled receptor (GPCR) that contains 11-ds-retinal as 
its intrinsic chromophore ligand, and it is highly concentrated 
in the ROS discs (Liang et al 2003; Nickell et al, 2007). Due 
mainly to its high endogenous expression, rhodopsin was the 
first structurally resolved mammalian GPCR (Palczewski et al, 
2000). In disc membranes, rhodopsin is tightly packed into 
paracrystalline dimer arrays, enabling optimal association 
with the heterotrimeric G protein transducin as well as with 
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additional regulatory components (Filipek etal, 2004; Fotiadis transducin, which in turn activates phosphodiesterase 6 
et al, 2004; Ciarkowski et al, 2005). Photon-activated (PDE6), leading to hydrolysis of cGMP and closure of the 
rhodopsin promotes the activation of the associated G protein cGMP-gated channels. This initiates ultra-fast phototransduc- 
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tion (Hamer et al, 2005), translating light energy first into a 
biochemical signal, followed by an electrical cue that is 
transmitted through the neuronal network of the retina. 
Adaptation to different light conditions, and regeneration of 
rhodopsin, is regulated at multiple levels, including through 
differential phosphorylation, differential calcium concentra- 
tions, and regulated enzymatic cycles, e.g., when regenerating 
11-ds-retinal (Lamb and Pugh, 2004). Disruption of these 
highly organized structures and processes by germline 
mutations can cause severe blinding diseases, such as retinitis 
pigmentosa (RP), rod-cone dystrophies, and congenital 
stationary night blindness (Berger et al, 2010) . 

Proteomic analyses of purified ROS have identified about 
500 proteins (Figure IB) that include metabolic enzymes, 
transport proteins, cytoskeleton elements, regulatory proteins, 
scaffolds, and housekeeping components, providing a detailed 
description of the outer segment protein repertoire (Kwok et al, 
2008; Figure 1C). In addition, the relative abundance has been 
determined for 150 proteins (Kwok et al, 2008) . Finally, studies 
have demonstrated that some of these proteins, such as 
arrestin, transducin, guanylate cyclase, and RhoA, localize 
differentially in the outer and inner segment in response to 
light/dark cycles (Hallet et al, 1996; Reidel et al, 2008; 
reviewed in Artemyev, 2008). To date, however, this wealth 
of data has not been fully analyzed nor integrated on a 
functional proteome-wide scale. It is anticipated that a large 
part of ROS proteins, and a core of the functional modules, will 
be common to all cells, while others will be photoreceptor 
specific (Hofmann et al, 2006) . Yet, similar to an assembly of 
music instruments, proteins organized as molecular machines 
can function in different, context-dependent ways. Connectiv- 
ity, as well as the timing and tuning of different modules, 
appears to be crucial for the proper orchestration of signal 
transduction as well as for parallel signal processing in a 
concerted manner: at the systems level, this results in the 
music of life. Although many details of the core phototrans- 
duction processes have been established, and mathematical 
models have been proposed (Dell'Orco et al, 2009), the overall 
orchestration of the outer segment functions, which include 
processes like disc shedding and renewal, protein transport 



along cilia, and light adaptation, are far from being under- 
stood. Further, how a variety of mutants of proteins primarily 
localized in the outer segments can cause visual impairment 
by perturbing the function of this organelle can only be 
speculated. For instance, a mutation could impair not only the 
proper folding of the protein but also its interactions with its 
partners within the physiologically functional protein net- 
works. Therefore, identifying protein interactions and their 
networks is an important step toward improving our under- 
standing of the molecular defects that underlie genetically 
inherited and age-related blinding diseases, and may directly 
lead to identifying novel disease-associated genes. 

There are several aspects that cannot be clearly determined 
through large-scale studies, such as the network dynamics, the 
simultaneous regulation of several distinct higher order 
biological outputs by one network, and the possibility that 
interactions detected for a particular protein might not be 
compatible simultaneously (Ito et al, 2001; Gavin et al, 2002, 
2006; Ho et al, 2002; Rual et al, 2005; Stelzl et al, 2005). As a 
consequence, information about the dynamics and the tempo- 
spatial resolution of networks has been limited to smaller 
signaling modules, such as receptor-initiated signal transduc- 
tion (Olsen et al, 2006; Becker et al, 2010). Structural 
information can help discriminate between direct and indirect 
interactions in a given complex. More importantly, it can add a 
dynamical value to the classical interaction networks by 
determining if two or more predicted partners of any given 
protein or complex can simultaneously bind to a target, or if 
they instead compete for the same interaction surface (Kim 
et al, 2006; reviewed in Campagna et al, 2008 and Kiel et al, 
2008). Integrating interaction data with protein expression 
information may assist in adding a dynamic dimension, and 
therefore a more realistic view, to the abstract 'organism 
interactome' (Hofmann et al, 2006) . 

Here, we have combined experimental data, literature 
mining, and structural information to provide a comprehen- 
sive view of the signal transduction network centered on 
rhodopsin (see the flowchart in Figure 2). Integrating 
structural information with the relative estimates of expres- 
sion levels allowed us to distinguish between mutually 



Figure 1 Proteomic description of the retina ROS inventory and GO analysis. (A) Schematic model of a rod photoreceptor cell (left) and its corresponding location 
within the retina (depicted in the micrograph to the right). Segments labeled in the model are ROS with enclosed stacks of discs membranes containing the visual pigment 
molecules rhodopsin; CC; RIS containing mitochondria, Golgi, and ER membranes, and vesicles in which opsin molecules are assembled before transported to the outer 
segment; and the cell body containing the nucleus and a synaptic termini, where neurotransmission to second-order neurons occurs. The micrograph depicts the vertical 
porcine retina with its cytoarchitectural organization labeled as photoreceptor outer segments (OSs); the outer nuclear layer (ONL) containing cell bodies of rods and 
cones; the outer plexiform layer (OPL); the inner nuclear layer (INL); the inner plexiform layer (IPL), and the ganglion cell layer (GCL). The retinal pigment epithelium 
(RPE) is localized above the photoreceptor cell layer (for details, see http://webvision.med.utah.edu). Retinal cells nuclei were stained with DAPI (magnification x 40). 
Insets show micrographs of the OS immunolabeled with anti-rhodopsin with an FITC-conjugated secondary antibody (magnification x 40; top inset) and of the OS 
preparation (magnification x 40; bottom inset). (B) Comparison of different proteomic data sets determined in ROS, based on proteins and the protein overlap identified 
in the proteomic analysis from this work and that of Kwok et al (2008). The union of the two data sets was defined as the initial experimental ROS proteome. 
(C) Functional modules and GO analyses of the filtered core ROS proteome. By performing an automatic and a manual GO search (based on the UniProt and KEGG 
databases), we characterized the 355 proteins (see Supplementary Table S2) to be involved in vision, signaling, transport, and channels (56), disc structure and 
morphology (7), housekeeping functions (73), cytoskeleton and polarity (67), vesicle, structure, and trafficking (60), and metabolism (92). Sub-modules/sub-functions of 
the GO terms are indicated as described in Supplementary Table S2 (1A, phototransduction/channels (33); 1B, retinol recycling (5); 1C, calcium signaling (18); 2A, disk 
morphology (2); 2B, link to ECM (5); 3A, protein folding (8); 3B, chaperones/heat shock (25); 3C, ubiquitination/degradation/proteasome (10); 3D, scaffolds/adaptor 
proteins (7); 3E, oxidative stress/cell redox homoestasis (9); 3F, apoptosis (2); 3G, others (2); 3H, signaling (10); 4A, regulation of cytoskeleton (34); 4B, cytoskeleton 
proteins (21); 4C, motor proteins (7); 4D, protein transport (1); 4E, axon guidance (4); 5A, endocytosis (10); 5B, exocytosis (8); 5C, Golgi endosome (11); 5D, vesicle 
transport/fusion (12); 5E, Golgi/ER/trafficking (19); 6A, glycolysis (20); 6B, tricarboxylic acid (5); 6C, ATP synthesis (25); 6D, lipid/fatty acids metabolism (9); 6E, amino- 
acid metabolism (9); 6F, one-carbon metabolism (4); 6G, nucleotide metabolism (6); 6H, glucose/lipid/phosphate/amino acid/ion transport (8); 61, pentose phosphate 
shunt (1); 6J, mevalonate (1); and 6K, others (4)). 
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compatible or mutually exclusive interactions, enabling us to 
structure a network of nodes and edges toward sub-networks 
and functional modules. The resulting network offers an 
unprecedented view of signal transduction in vision and 



suggests a light-dependent orchestration of the core vision 
pathway to functions that have so far not been related to this 
pathway, such as cytoskeleton dynamics, vesicle transport, 
and energy metabolism. The specific light-dependent con- 
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Figure 2 Experimental and computational workflow. The flow charts of experimental (yellow boxes) and bioinformatic (green boxes) methods used in this work are 
shown. The initial ROS proteome was generated based on the union of proteins identified in bovine ROS in this work and those from a proteomic analysis of porcine ROS 
(Kwok et al, 2008). After filtering, a high-confidence ROS proteome was defined. A static ROS interactome was compiled by literature mining. In addition, new 
experiments were performed in ROS in this work (co-sedimentation and co-IP). Further, we performed structural analyses and homology modeling to distinguish 
between compatible and mutually exclusive interactions. This enabled us to break the network of nodes and edges into functional machines or sub-networks and 
modules. The comprehensive multiscale network highlights new predicted links and functions. Finally, disease-associated genes were identified and modeled into 
available structures. 
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nectivity of rhodopsin to these functions is likely to be 
conferred by small GTPases and their regulators and interact- 
ing proteins, such as the prenyl-binding protein PDE5 and 
other prenylated proteins. This would establish a dynamic and 
light-dependent mode of regulating the localization (to the 
membrane or cytosol), and thus the activity, of these GTPases. 

Results 

ROS proteome determination and contaminant 
removal 

To determine the proteomic content of photoreceptor outer 
segments, dark-adapted porcine ROS and outer segments discs 
were isolated as previously described (Swiatek-de Lange et al, 
2008; see Materials and methods) . Proteins were then resolved 
by one-dimensional gel electrophoresis (IDE) and identified 
by mass spectrometry (MS). In three independent experi- 
ments, a total of 50 proteins were identified by MALDI-TOF 
from both ROS and ROS discs, and 434 proteins by the more 
sensitive Orbitrap-LC-MS/MS (Supplementary Table SI). The 
union of the two data sets resulted in a total of 444 proteins, of 
which 410 could be mapped to the human proteome. 
Comparing our data set to a recent proteomic study that 
identified 516 proteins from bovine ROS (of which 487 mapped 
to the human proteome; Kwok et al, 2008) resulted in an 
overlap of 217 proteins (Figure IB). We then created a unified 
data set consisting of 680 human proteins, defined as our 
'initial experimental ROS proteome'. 

To further refine the protein list presented in the initial 
experimental ROS proteome, we applied heuristic filtering 
procedures to remove proteins that might have contaminated 
the ROS fraction in either experiment (from the surrounding 
cells or from other cellular domains of the rod photoreceptor). 
First, we looked at the functional annotations of the initial 
experimental ROS proteome to identify annotations that 
contrasted with the expected properties of an ROS protein 
(such as transcription factors, nuclear proteins, and mitochon- 
drial proteins), based on GO terms from the UniProt database 
(Supplementary Table S2). Second, we performed a detailed 
manual functional analysis based on UniProt, the KEGG 
database, and relevant literature. This revealed 81 putative 
contaminants (Supplementary Table S2); of these, 68 were 
found in only one of the two experimental data sets (e.g., ours 
and that from Kwok et al, 2008), further supporting their 
classification as contaminants. The 13 proteins identified in 
both sets are synaptic proteins and G proteins believed to be 
expressed only in cones (Kwok et al, 2008) . We, thus, retained 
605 proteins after these analyses. 

We next removed all proteins for which there was no 
interaction data or further experimental evidence about their 
presence in ROS (protein group 1; Supplementary Table S2); 
this information was compiled from the literature and 
databases such as MINT (for details, see Materials and 
methods and http://mint.bio.uniroma2.it/mint/) (Zanzoni 
et al, 2002; Chatr-aryamontri et al, 2007; Ceol et al, 2010; 
Supplementary Table S3). A total of 347 proteins passed this 
filter. From the analysis of published information, eight 
additional proteins that lacked associated protein-protein 
interactions (PPIs) were nevertheless considered to be bona 

© 2011 EMB0 and Macmillan Publishers Limited 



fide ROS proteins with important functional and/or structural 
roles and were thus retained in the ROS proteome; these were 
the retinal-specific ATP-binding cassette transporter ABCA4, 
the cellular retinaldehyde-binding protein RLBP1, the photo- 
receptor outer segment alRrans-retinol dehydrogenase RDH8, 
peripherin-2 (PRPH2), the ROS membrane protein ROM1, 
RabllB, RP1, and fascin 2 (FSCN2). This filtering procedure 
left us with 355 bona fide ROS proteins or proteins that are 
dynamically localized to ROS (protein groups 2 and 3, 
respectively, in Supplementary Table S2). These proteins 
represent the 'core ROS proteome'. 

Functional modules of the core ROS proteome 

We next classified the core ROS proteome into six functional 
groups based on the above information, and we annotated 
lipid modifications, such as prenylation and geranylation 
(Supplementary Table S2; Figure 1C): (1) vision, signaling, 
transporters, and channels: 56 proteins have functions that are 
either directly associated with vision or support visional 
functionality (i.e., visual cycle, protein homeostasis, or energy 
production). This module contains well-known members of 
the phototransduction pathway, including the core signal 
transduction of light (Dell'Orco et al, 2009) and the visual cycle 
involved in regenerating ll-d5-retinal necessary to comple- 
ment photosensitive rhodopsin after photo-bleaching (Lamb 
and Pugh, 2004) . We further included here proteins involved 
in Ca 2 + -dependent signaling and proteins associated with ion 
channels that regulate photoreceptor membrane conductance 
and polarity; (2) outer segment structure and morphogenesis: 
the seven proteins in this group are those implicated in outer 
segment structure and disc morphogenesis (Molday et al, 
1987; Poetsch et al, 2001), and those that link the cytoskeleton 
to the extracellular matrix (ECM), such as oc and [3 catenin; (3) 
housekeeping: in this group of 73 proteins, we consider 
protein-folding chaperones and heat-shock proteins, members 
of the ubiquitination/degradation-proteasome machineries, 
scaffold proteins such as the 14-3-3 family members, and 
proteins involved in oxidative stress, cell redox homeostasis, 
and apoptosis regulation (De La Paz and Anderson, 1992; for 
review, see Wenzel et al, 2005); (4) cytoskeleton and polarity 
(67 proteins): this group contains cytoskeleton proteins, such 
as actin and tubulin, as well as their respective binding 
proteins and molecular motors, proteins involved in regulating 
cytoskeleton dynamics including GTPases, and intermediate 
filaments. Many of these are associated with the CC and the 
axoneme of ROS and might therefore be present at low 
concentrations (reviewed in Adams et al, 2008). Proteins that 
are known to function in axon guidance were also added to this 
class; (5) vesicles formation and trafficking: we included here 
60 proteins involved in Golgi function, protein and vesicle 
transport, and fusion, as well as the annexins that function in 
exocytosis and phagocytosis; (6) metabolism: we included 
here 92 proteins related to metabolism, in processes such as 
glycolysis, ATP synthesis, nucleotide, and fatty acid and 
carbon metabolism. Interestingly, we found that several of 
these are metabolic proteins involved in energy production 
(about 50 % of the enzymes detected in this group are involved 
in glycolysis and about 20% in the tricarboxylic acid path- 
way), including ATP synthase, the activity of which has 
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recently been demonstrated in intact discs (Panfoli et al, 2009) . 
This suggests that the ATP used in vision signaling is indeed 
produced within ROS and is probably fueled by glucose 
transported along the cilium. Indeed, this study identified a 
glucose transport protein, SLC2A1, as an ROS protein, 
supporting this hypothesis. High energy demands of ROS, 
and a capacity of only limited diffusion through the 
interconnecting cilium, may require on-site production of ATP. 



Network reconstruction and structural modeling of 
the ROS interactome 

Information about the PPIs among the core ROS proteome was 
mined from protein interaction databases to assemble an ROS 
protein network (Supplementary Table S3; Zanzoni et al, 2002; 
Chatr-aryamontri et al, 2007; Kerrien et al, 2007; Ceol et al, 
2010). The protein interaction degree ranges from 1 to 179, 
with the highest number of interaction partners for actin, 
tubulin, 14-3-3 family members, heat-shock protein members, 
and ERK (Supplementary Figure SI A). 

Overall, the complete core ROS proteome PPI network 
consists of 5337 interactions among its members (Supplemen- 
tary Table S3). The experimental evidence for most of these 
interactions (5047) was based on co-immunoprecipitation (co- 
IP) or pull-down experiments, which offers little support for 
their direct nature (Supplementary Figure SIB). In addition, 
many of the edges in the network are supported by single 
experimental pieces of evidence (>85% of the PPI), often 
derived from high-throughput approaches. Thus, we refer to 
this network, which represents all the interactions that we 
could retrieve from published data, as a 'fuzzy ROS inter- 
actome\ since it contains many interactions supported by only 
one non-binary piece of evidence. 

Next, we aimed at increasing the information content of the 
network by structural modeling. Pairs of interacting proteins 
often share common structural features with other interacting 
pairs of known structure (domains and linear motifs) . We use 
structural information, combined with computational tools, to 
support low-confidence experimental interaction evidence and 
to determine whether two interactions involving a common 
partner are compatible or mutually exclusive (see Supplemen- 
tary Material 1 and Supplementary Figures S2 and S3). We 
considered two levels of structural evidence that may support 
any given interaction. First, for each pair of members of the 
core ROS proteome, we searched the PDB database (http:// 
www.pdb.org) for protein complexes of known structures 
whose elements share at least 70% homology with the query 
proteins. By this approach, we identified 84 complexes in the 
ROS core network whose structures could be confidently 
modeled on homologous structures (Supplementary Table S4) . 
Most of the interactions for which there are X-ray structures, or 
structures from close homologs, are found between the 
connecting proteins in the modules 1 or 4 [vision, signaling, 
transporters, and channels and cytoskeleton and polarity) as 
well as among proteins involved in interactions connecting 
these two modules (Supplementary Figure S4) . 

Next, using a lower level of structural detail and confidence, 
we exploited the notion that similar domain pairs are likely to 
interact in a similar way ('nature repeats itself) (Aloy and 
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Russell, 2002). For example, members of the Ras family and 
proteins containing a Ras-binding domain (RBD) are likely to 
use the same interaction surface when they interact (see Kiel 
and Serrano, 2006). To overlay a domain-level model on the 
ROS network, we represented each of the 355 nodes as a stack 
of Pfam domains (http://pfam.janelia.org/; Supplementary 
Table S5). We then searched the 3Did database (http:// 
3did.irbbarcelona.org/; Stein et al, 2005, 2009) for structural 
evidence of pair-wise interactions between any of the domains 
in our database (for details, see Supplementary Material 1). 
Structural evidence was found for 352 pair- wise interactions, 
excluding pairs that had already been identified by comparison 
with homologous crystallized complexes (Supplementary 
Table S4). A confidence 'interaction score' of ^2.3 for the 
identification of the interacting pairs was obtained by 
interrogating the InterPReTS server (http://www.russelllab. 
org/cgi-bin/tools/interprets.pl/; Aloy and Russell, 2003). This 
score was validated using a yeast two-hybrid positive and 
negative binding data set described by Vidal and co-workers 
(Rual et al, 2005) . We found a confidence of over 70 % that two 
proteins containing the target domains will interact in a two- 
hybrid experiment when the InterPReTS score was >2.3 (see 
Supplementary Material 1, Supplementary Figures S2 and S3). 
Of the 352 interactions that had a hit in the 3 Did database, 107 
had InterPReTS scores higher than our chosen threshold and 
were therefore annotated as 'supported by structural evi- 
dence'. A total of 191 interactions supported by structural 
evidence (that is, the 84 interactions with known or closely 
related structures, and the 107 with significant InterPReTS 
scores) were merged with the literature-based interaction 
network. Interactions that could be annotated with structural 
evidence were mainly found within the functional modules 
(Supplementary Figure S4) . 

To increase the confidence in the resulting network, edges 
that were only supported by a single piece of evidence from 
any type of experiment except yeast two-hybrid experiments 
were removed (Supplementary Table S6) , with the exception of 
interactions for which there was also structural information 
available (i.e., a three-dimensional structure of the complex 
itself or of a highly homologous complex) . This curated static 
network {'high-confidence ROS interactome') comprises 660 
edges and links the majority of the nodes (with 266 proteins, as 
indicated in Supplementary Table S2 and Figure 3A) that were 
present in the original network. The missing nodes are equally 
distributed among the proteins with respect to their GO terms, 
although an enrichment for proteins assigned to the classes 
retinol recycling and metabolism was observed (of 80 and 50 % , 
respectively) . 

By considering only edges supported by at least one 
evidence of direct binary interaction, we obtained a 'high- 
confidence binary ROS interactome' that contains 222 nodes 
(note that most of the nodes that were not captured by this 
network are annotated with metabolism ontology terms), 
linked by 349 edges (indicated as binary in Supplementary 
Table S6). Except for reactions involving guanylate kinase 
(GK) and nucleoside diphosphate kinase (NDPK), a nucleotide 
and cyclic nucleotide modifying enzyme, all interactions of the 
core vision pathway (Ridge et al, 2003; Wensel, 2008; 
Dell'Orco et al, 2009) are represented in our PPI network as 
true binary interactions. Thirty-five nodes have more than ten 
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Figure 3 The high-confidence ROS interactome and the high-confidence binary ROS interactome. (A) The high-confidence ROS interactome. The 660 higher 
confidence interactions of the ROS interactome are listed (Supplementary Table S6). The size of the nodes indicates the number of interaction partners for a given 
protein (of > 10 or >20). Edges with binary evidence are indicated with blue, while edges supported by more than one piece of evidence are indicated in gray. Proteins 
are colored according to their function. (B) The high-confidence binary ROS interactome. Modules and sub-modules are shown, and only the interactions of proteins 
from two different modules are indicated (see Supplementary Material 2). The number of proteins implicated in diseases in each category is indicated. 
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interaction partners, with a maximum degree of 55. Only 2 
proteins in the vision category have > 10 interaction partners 
(CALM1 and CAMK2 A) . Most of the interactions involve the 
heat-shock proteins, 14-3-3 family members, ERK (MAPK1), 
tubulins, and actin. More than ten interaction partners were 
found in the metabolism branch for glyceraldehyde-3 -phos- 
phate dehydrogenase (GADPH) and for the two ATPase 
subunits. 

Additionally, 109 direct binary interactions connected the 
defined functional modules, and 240 binary interactions 
connected proteins within modules (Figure 3B; for a detailed 
description of the connections between the modules, see 
Supplementary Material 2). Out of those 240, 188 are classified 
within sub-boxes/sub-functions. The observation that roughly 
two thirds of the interactions were found within functional 
modules, and only a third between modules, provides 
confidence to our module classification and manual functional 
annotation. Interestingly, the most highly connected modules 
are module 1 (vision) and module 4 (cy to skeleton), illustrating 
the important crosstalk between the core vision pathway and 
the cytoskeleton. The less connected modules are the ones 
involved in the structure of the discs and in metabolism. As 
expected, the housekeeping module, despite having fewer 
connections, is linked to all other modules. 

The high-confidence ROS interactome suggests 
new functional links 

Using our curated high-confidence binary ROS interactome as 
a basis, we decided to analyze in more depth the core vision 
pathway, which is probably one of the best-studied biochem- 
ical pathways (Ridge et al, 2003; Wensel, 2008; Dell'Orco etal, 
2009). We extended the published core vision pathway 
(Dell'Orco et al, 2009) using evidence from our high- 
confidence network and indicated structural coverage and 
outputs to different functional cellular processes emanating 
from the proteins in the pathway (Figure 4; for a detailed 
description, see Supplementary Material 2). Of these, we 
decided to validate the link to the GTPases RhoA and Racl 
(Figure 4, link A and link D), which suggests a link between 
vision activation and cytoskeleton reorganization. 

Rho-Rac1 and the cytoskeleton connection 

Previous work has demonstrated functional links between 
rhodopsin, certain GTPases (Mitchell et al, 1998) (most 
prominently transducin), and the cytoskeleton. S-arrestin 
specifically binds to activated and phosphorylated rhodopsin, 
inhibiting activation of transducin and terminating photo- 
transduction (Kiihn, 1978; Kiihn et al, 1984; Wilden et al, 
1986). Nair et al (2004) have shown interactions between 
S-arrestin and microtubules (Figure 4, link G) . 

We were able to confirm that small GTPases Rac and the 
GTP-bound form of RhoA bind rhodopsin, as has been 
previously described (Wieland etal, 1990a, b; Balasubramanian 
and Slepak, 2003; Gray et al, 2008; Figure 4, link A). For this, 
we performed co-segregation/co-sedimentation experiments 
to reveal proteins within large complexes, as described 
previously to analyze the light-harvesting complex of photo- 
system II in plants (Swiatek-de Lange et al, 2008; Supplemen- 
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tary Figures S5 and S6A). These experiments indicated that 
Racl, Rho, and CRMP-2 were present in a large complex that 
also contained cytoskeletal proteins, rhodopsin, and compo- 
nents of the vision pathway. Although the low resolution of the 
technique, and the complexity of the patterns, preclude using 
these experiments to add new binary interactions to the ROS 
network, it can be used to corroborate interactions supported 
by further experimental evidence or from the literature 
(Supplementary Tables S3 and S6). Using BN-PAGE (Schagger 
and von Jagow, 1991; Nijtmans etal, 2002; Camacho-Carvajal 
et al, 2004) or IP experiments in combination with either MS or 
subsequent immunoblotting, we obtained further evidence for 
the existence of large complexes containing rhodopsin that 
also included the cytoskeletal proteins actin and tubulin as 
well as its specific regulators such a RhoA, Racl, and CRMP-2 
(Figure 5; Supplementary Figure S6B). 

We next isolated the protein partners of the glycosylated 
N-terminus of rhodopsin by using concanavalin A affinity 
purification (Plantner and Kean, 1976; De Grip, 1982). The 
interactions between rhodopsin, RhoA, and CRMP-2 were 
confirmed by these concanavalin A pull-down experiments, 
and in part by additional co-IP experiments in which we 
detected rhodopsin associated with the core signaling com- 
plexes of the visual pathway including transducin and, again, 
with Rho, Racl, and CRMP-2 (Supplementary Figure S6C and 
D). To confirm that CRMP-2, Racl, and ROCK II are indeed 
bona fide ROS proteins rather than contaminants, we 
performed immunohistochemistry for these proteins. Indeed, 
all three proteins were constituents of ROS on cryosections of 
porcine retina (Figure 6) . Despite considerable efforts, we were 
not able to confirm the presence of RhoA due to a lack of 
selectivity of various antibodies against RhoA in retinal 
sections. 



Functional analysis of the PDE5-Rac1 complex 

PDE5 has been reported to bind to prenyl-modified proteins, 
such as several small GTPases and rhodopsin kinase (Hanzal- 
Bayer et al, 2002; Zhang et al, 2004), and it appears as an 
important node within our network (Figure 4, link D). PDE<5 
could thus have a critical regulatory role both in facilitating the 
transport of prenylated target proteins along the cilia together 
with Arl3 (Figure 4, link C; Veltel and Wittinghofer, 2009) and 
in serving as an effector or guanine nucleotide dissociation 
inhibitor (GDI) for many GTPases, such as Arf, Racl, RhoA, 
and Rab, all of which are expressed in photoreceptors 
(Figure 4, link D). Therefore, we tested whether PDE<5 
functions as a GDI for Racl in ROS. First, we demonstrated 
that PDE5 and Racl are colocalized in ROS using immunohis- 
tochemistry (Figure 6; Supplementary Figure S7). Second, we 
showed that PDE5 and Racl colocalize in ROS in native protein 
complexes, by using dark-adapted ROS separated by BN-PAGE 
(Figure 7A). In dark-adapted ROS, PDE5 was part of distinct 
complexes that ranged from high molecular weight complexes 
of 660 kDa to smaller complexes of around 90kDa; and 
interestingly, Racl colocalized with PDE<5 within different 
complexes between the soluble and membranous fractions. 
Third, we tested whether PDE5 could dissociate Racl from 
ROS membranes in vitro (see Materials and methods) . Indeed, 
adding recombinant human (rh) PDE5 led to the solubilization 
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of Racl from the ROS membranes (Figure 7B) . Solubilization 
occurred in a dose-dependent manner with increasing 
amounts of rhPDE<5. As a positive control, we verified that 
PDE5 solubilized PDE(3 from ROS membranes, as previously 
described (Florio et al, 1996). Thus, PDE5 can solubilize Racl 
from ROS membranes, a feature characteristic of GDIs. 

All of the interactions determined here — with the exception 
of the ones identified by co-sedimentation, as this method is 
considered as weak evidence for physical interactions — were 
added as supporting evidence to our network (Supplementary 
Table S6) . In total, co-purification and co-elution experiments 
supported 60 interactions that had been included in our 
network based on the literature, and new evidence for 175 
interactions from our co-IP results was added. Additionally, 
our results supported five interactions that had structural 
evidence (with INTERPRETS score ^2.3). 

Restricted to a single new pathway (Racl/RhoA- 
PDE5— CRMP-2), our experimental data support the physiolo- 



gical relevance of our network. It should be noted, however, 
that these data cannot be considered to be complete or free of 
false positives, since the number of interactions tested and 
validated was small when considering the extent and complex- 
ity of the network. 



Discussion 

In this work, we investigated the protein interaction network 
in a highly specialized cellular region of the mammalian 
photoreceptors, the ROS. Graphs representing protein interac- 
tions are idealized descriptions of all the interactions that can 
possibly occur in an organism. The realization that, in any 
given cell type, only a fraction of these interactions can 
possibly occur prompted the development of approaches to 
combine different genome-wide information to build interac- 
tion networks that are either specific for a cell type (Bossi and 
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Figure 5 Graphical representation of experiments performed in this work and its comparison with interactions described in the literature. Protein complexes that were 
obtained using Rac1 , RhoA, or Rac1 as the bait protein are displayed within orange, blue, and yellow circles, respectively. The Rac1 and RhoA complexes were identified 
by western blot, and the Rac1 complex by Orbitrap. The overlap of the three circles indicates the proteins that were identified in the same complex in one of the three 
experiments. Connecting lines between proteins indicate either binary or co-IP interactions from the literature, or from BN-PAGE or co-sedimentation interactions as 
determined in this work. Proteins are colored according to their function. 



Lehner, 2009) or that change dynamically, such as during the 
cell cycle or after specific pathways have been induced. Here, 
we take this one step further and propose a protein interaction 
network for a structurally very distinct and functionally highly 
specialized region of the mammalian photoreceptors: the ROS. 
In addition to proposing novel interactions, we present a 
structural model that allows us to discriminate between 
protein interactions that are compatible and those that are 
mutually exclusive. 



A curated and structure-based PPI network central 
to rhodopsin 

We first generated a ROS-specific protein interaction network 
by combining proteomic expression levels in ROS with 
interaction information, which we mined from the literature 
and then subsequently supplemented with our new data 
pertinent to the description of protein complexes in their 
physiological context. We next performed structural analysis 
of the curated network by decomposing proteins within this 
network into domains. This step allowed us to validate 
interactions at a domain level and to thereby increase the 
confidence in the network. By reassembling the decomposed 
network based on structural constraints into structure- 
functional modules, we were able to define logical relation- 
ships between the network nodes, and to define sub-networks 
that physically and functionally fit into molecular machines. 
Finally, we annotated these functional modules according to 
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their respective physiological processes to derive a network of 
pathways and processes. Based on a compilation of experi- 
mental evidence and several layers of expert as well as 
automated curation, filtering, and modeling, the resulting 
network represents a multiscale description of wiring and 
physical connectivity in the ROS of photoreceptors. The 
extended core pathway shows how rhodopsin activation- 
deactivation leads to other possible functional effects in 
addition to its primary function of signaling for closing the 
cGMP-gated cation channel. Thus, in addition to its relation- 
ship with the module of (1) vision, signaling, transporters, and 
channels, wiring rhodopsin to (2) outer segment structure and 
morphogenesis, (3) housekeeping, and (4) cytoskeleton and 
polarity suggests a regulation of cytoskeleton assembly-dis- 
assembly and dynamics, vesicle and Golgi trafficking, and 
transport along the interconnecting cilium of photoreceptors 
by rhodopsin. Connections between active rhodopsin and Arf4 
(Deretic et al, 2005; Mazelova et al, 2009), and between PDE5 
and Rabl3 and the GTP-bound form of Arl3 (Hanzal-Bayer 
et al, 2002) , also link the vision cycle to vesicle trafficking and 
structure (Figure 4B and C) . We experimentally validated two 
of the proposed new functional links. Our results suggest a link 
between rhodopsin, Racl, RhoA, ROCK II, and CRMP-2. This 
points to a second, not yet experimentally tackled pathway 
that is influenced by light, which appears to be a delineation of 
an archetypical G protein-regulated pathway known to be 
active in growth cone dynamics and collapse (Liu and 
Strittmatter, 2001). RhoA binds to CRMP-2 (gene name 
DPYSL2, Figure 5), a scaffold protein involved in actin 
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Figure 6 Immunohistochemical analyses of porcine retina. Cryostat sections of the retina were stained with primary antibodies (red) against indicated proteins, and 
nuclei were counterstained (blue). The images on the left were taken from the outer retina (outer segments (OS), inner segments (IS), outer nuclear layer (ONL), and 
outer plexiform layer (OPL)). Images in the middle are an overlay of antibody staining, nuclei staining, and DIC optics (Nomarski). Images on the right were taken with 
higher magnification, to focus on the OS and IS. All indicated proteins were unambigiously identified as constituents of ROS. Control sections without primary antibodies 
showed no staining (Supplementary Figure S7). 

cytoskeleton dynamics in neurons that regulates growth cone axon guidance, where it dynamically regulates the antagonis- 
dynamics. CRMP-2, working through the GPCR lysophospa- tic effects of RhoA and Racl. Regulated by a Rho-associated 
tidic acid receptor, has been described as a crucial molecule in kinase (ROCK), CRMP-2 promotes either outgrowth or 
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Figure 7 Experimental evidence that PDE5 acts as a GDI for Rac1 in ROS. 
(A) PDE5 and Rac1 colocalize in ROS in native protein complexes. After 
solubilization with p-DM, native ROS protein complexes from soluble and 
membranous fractions of light- and dark-adapted ROS were separated by 
BN-PAGE. Components of the native protein complexes were separated by 
SDS-PAGE for second-dimension electrophoresis. Western blots with anti-Rad 
and anti-PDE5 antibodies showed that PDE5 and Rac1 colocalized but were in 
different complexes in ROS depending on the dark-adapted state of the retina. 
Colocalization of PDE5 and Rac1 seemed to be stronger in the dark-adapted 
state, where both proteins colocalized to the soluble and membranous fractions. 
In light-adapted ROS, colocalization of PDE5 and Rac1 was detected only in the 
membranous fraction but not in the soluble fraction. (B) In-vitro solubilization of 
Rac1 GTPase from light- and dark-adapted ROS membranes. Membranes 
isolated from light- or dark-adapted ROS were incubated for 1 h at 37°C with 
different amounts of recombinant human PDE5 (rhPDES) or buffer alone, and 
the unsolubilized material was recovered by ultracentrifugation. Immunoblots 
with anti-Rad or anti-PDES antibodies showed that PDE5 solubilizes Rac1 from 
ROS membranes in a dose-dependent manner. Since it has been previously 
determined that PDE5 solubilizes PDE5 from ROS membranes in a dose- 
dependent manner (Florio et al, 1996), this was used here to demonstrate the 
functional activity of the rhPDES protein. 



collapse in response to active RhoA or Racl, respectively (Hall 
etal, 2001). When RhoA-GTP levels are high, more CRMP-2 is 
phosphorylated by the Rho-effector kinase ROCK, and thus 
less non-phosphorylated CRMP-2 is complexed with Racl, 
leading to cytoskeleton collapse (reviewed in Liu and 
Strittmatter, 2001). CRMP-2 can bind directly to tubulin 
heterodimers to promote microtubule assembly (Fukata 
et al, 2002) . This presents the exciting possibility that GPCR 
rhodopsin autoregulates its own axonal/dendritic guidance 
and possibly regulates outer segment growth via the arche- 
typical mechanisms of axon guidance. Based on this scenario, 
the outer segment would function as a continuously extending 
growth cone, autoregulated by light and other as-yet uni- 
dentified guidance cues that may be produced in other retinal 
cells, most notably in the retinal pigment epithelium. 

We additionally provide experimental evidence that PDE5 
could act as a GDI for the small GTPase Racl. PDE5 could thus 
have a very crucial regulatory role: (1) in transporting 
prenylated target proteins (Zhang et al, 2004) along the cilia, 
together with Arl3 (Veltel and Wittinghofer, 2009) and (2) as 



an effector or GDI for many GTPases (Hanzal-Bayer et al, 
2002), such as Arf, Racl, RhoA, and Rab, by keeping them GDP 
bound and inactive. This is important since we did not find the 
conventional RhoGDI in ROS, suggesting that PDE5 could 
indeed substitute for this function in ROS (similar to that 
demonstrated for the small GTPase Rab 13 in ROS; Marzesco 
et al, 1998). However, despite the structural similarity of the 
PDE5 and RhoGDI domains (Scheffzek et al, 2000; Hanzal- 
Bayer et al, 2002), we learned by superimposing the Racl- 
RhoGDI with the Arl2-PDE<5 structure that these two interac- 
tions depend on different moieties for binding (Supplementary 
Figure S8) . We provide experimental evidence in this work that 
PDE5 could act as a GDI for Racl. We did not find any GEFs or 
GAPs for small GTPases in our network but only GDIs 
(ARHGDI for RhoA, PDE5 for Racl, and GDI1 and GDI2 for 
Rab proteins). Interestingly, this could suggest that these are 
not regulated by the usual switch-like mechanism of GTPase 
regulation, but rather by a gradient activation, in which the 
activity of active RhoA is determined only by the concentration 
of RhoGDI, keeping RhoA in the inactive form. 



The role of Ca 2+ in vision cycle, phototransduction, 
and actin cytoskeleton changes 

Intracellular Ca 2+ concentrations influence the activities of 
numerous kinases, such as different PRC isoforms, the PKA 
kinase, Ca 2 + /calmodulin-dependent kinases, and the two 
CaMK-II isoforms, all of which are integral to the network. 
Predicted kinase phosphorylation sites from CaMK-II, PKA, 
PKC, MAPK, and PKD are summarized in Supplementary 
Table S7. Several Ca 2 + -regulated kinases phoshorylate cytos- 
keletal target proteins, such as actinin and myosins, and small 
GTPases and their regulators. This opens the intriguing 
possibility that the nucleotide state and the dynamic spatial 
cellular distribution of several small GTPases are controlled by 
Ca 2 + . As perturbed Ca 2 + homeostasis is a consequence of the 
activity of a perturbed visual pathway in specific forms of RP 
(Paquet-Durand et al, 2010), this is likely to affect a variety of 
critical pathways and thus generate a systemic perturbance of 
ROS physiology. Our network reveals several direct binary 
connections between Ca 2 + -regulated proteins and cytoskele- 
ton proteins: CaMK2A with actinin, calmodulin with GAP43 
(neuromodulin) and S1008 (tubulin polymerization initia- 
tion), and PKC with 14-3-3 family members. Calmodulin is 
known to have a wide range of effector binding specificity, 
which dynamically changes with Ca 2+ binding. Calmodulins 
1 and 3 were linked to about 10 proteins from the two modules, 
cytoskeleton and vesicle transport. Calmodulin (CALM3 or 
CALM1) can bind to the cytoskeleton regulator spectrin a, 
actinin (ACTN2 and ACTN4), and the myosin motor protein 
MY06. Therefore, calmodulin proteins could provide an 
important link between Ca 2 + signaling and regulation of the 
actin cytoskeleton, with spectrin having a critical role in 
organizing and maintaining membrane sub-domains that 
harbor rhodopsin (Berghs et al, 2000). Further, a Ca 2 + - 
dependent kinase, CaMK2A, was found to directly contact 
actinin-1, -2, and -4, and to be in a ternary complex with 
densin, a synaptic adhesion molecule (Walikonis et al, 2001), 
which is not present in our network (as it was not taken into 
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consideration). Another link appears between calmodulin and 
RalA and RalB, both of which are involved in trafficking: RalA 
has a role in exocytosis regulating exocyst assembly, while 
RalB interacts with EXOC8 (a part of the exocyst complex); 
RALBP1 is an effector of both RalA and RalB. Ca 2+ activity is 
also likely to regulate metabolic activities through IHD3A, 
recoverin, and neurocalcin: the hippo calcin-like protein 1 is a 
recoverin-like protein that was suggested to have an anti- 
apoptotic function and might protect photoreceptors from 
Ca 2 + -induced cell death (Krishnan et al, 2009). 

Finally, Ca 2 + could have an important role in the light-dark 
cycle by affecting PKA activity. Phosphorylation of RGS9-1 by 
PKA (Balasubramanian et al, 2001) is regulated by light and 
Ca 2 + , and results in the reduction of RGS9-1 GAP activity: with 
light, RGS9-1 causes rapid Toc-GTP inactivation and photo- 
receptor recovery, while in the dark, PKA is activated by rising 
concentrations of Ca 2+ and cAMP, which in turn phosphor- 
ylates RGS9-1. In this way, GAP activity is reduced, the active 
transducin lifetime is prolonged, and the photoresponse is 
strengthened (Balasubramanian et al, 2001). 

While it remains to be seen how all of these connections are 
orchestrated, and to which degree they impact vision home- 
ostasis, there is no doubt that Ca 2+ has a crucial role in ROS 
functionality. 



Structural information, structural coverage, and 
AND' and XOR' gates 

Structural information allows the confidence of any indepen- 
dent interaction evidence to be tested and at the same time can 
add topological information to the molecular level by defining 
sites or interaction domains. When several proteins can bind 
to a single protein, the various interactions can occur 
simultaneously or can be mutually exclusive (reviewed in 
Santonico et al, 2005 and Kim et al, 2006). If two or more 
proteins compete for the same binding site, then it seems 
unlikely that binding can occur simultaneously, whereas 
binding to topologically distinct sites may occur at the same 
time. At the level of graph representations within a network, 
structural information can thus support logical constraints. 
Here, it is important to mention that interactions were defined 
as exclusive or compatible from a structural point of view, and 
that this cannot be directly translated to biological terms in all 
cases (i.e., for competition to occur, the target protein should 
be at lower concentration than the competing ones). When 
both competitors are present at the same place and time, 
changes in concentration levels or additional regulatory 
constraints (e.g., those introduced by post-translational 
modifications) could regulate competition. 

Structurally superimposing domains onto interactions 
allowed us to define ternary complex formation and, 
importantly, to model both the composition of macromole- 
cular assemblies and its dynamic dissection into mutually 
exclusive complexes (Supplementary Figure S9). With this 
information, we can add dynamics to the network, using the 
following AND' and 'XOR' ('XOR'=exclusive OR) logical gate 
symbols: if three or more proteins can interact at the same 
time, then they are compatible (indicated with AND'), while if 
three or more proteins cannot interact simultaneously, then 
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they are mutually exclusive (indicated with 'XOR') (Figure 8). 
Competitors are frequently found in highly dynamic processes 
or may dynamically connect a given protein to different 
signaling and functional modules. The structural and interac- 
tion analyses of the core vision pathway and its cytoskeleton 
branch show several examples of non-compatible ('XOR') 
interactions (Figure 8) . For example, rhodopsin may interact 
with transducin, arrestin, or rhodopsin kinase (in the core 
vision pathway). It may also interact with Racl or RhoA 
(which are antagonists in cytoskeletal dynamics) or with Arf4 
(involved in trafficking). Changes in rhodopsin activation, 
concentration, and localization, or in its activation states, may 
therefore switch signaling into different pathways. Further, 
rhodopsin localization during ciliary transport and disk 
formation, and dynamic changes in concentrations and 
activation states in response to light, can alter the array of 
rhodopsin binding partners, since these are determined by the 
phosphorylation state of rhodopsin on the one hand and the 
availability or concentration of binding proteins on the other 
hand. Interestingly, AND' gates are mainly found in the 
housekeeping, structure and polarity, and metabolism 
branches, e.g., within large functional complexes, such as 
the T-complex, the proteasome, tubulin, and the ATP synthase 
machinery. 'XOR' gates, which are mainly prevalent in the 
vesicle structure and trafficking branch, indicate switch 
behavior or redundant protein functions, such as for Rab 
GTPases (Del Conte-Zerial et al, 2008). In the vision branch, 
both AND' and 'XOR' gates synergize. This may allow 
dynamic tuning of light and dark states. However, all 
connections from the vision module to other modules are 
'XOR' connections, suggesting that competition, together with 
local protein concentration changes, could be important for 
transmitting signals from the core vision module. 

The vision network and disease 

A large fraction of retinopathies involve the degeneration 
of rod photoreceptors; these include RP, syndromes incorpor- 
ating retinal degeneration with different associated pheno- 
types (such as Usher's syndrome and Bardet-Biedl syndrome), 
and Leber congenital amaurosis (LCA), a congenital form of 
retinal degeneration. An increasing number of genes and 
proteins has been implicated in these pathologies (http:// 
www.sph.uth.tmc.edu/Retnet/). These proteins include: 
components of the visual transduction cycle; structural 
components of the cytoskeleton, rod and/or cone photore- 
ceptor outer segment disc membranes; components of synth- 
esis and recycling of the retinoid; transcription factors 
(including CRX and NRL) and splicing factors; those involved 
in signaling and cilium maintenance, phagocytosis of the outer 
segment discs of the photoreceptors, and trafficking of 
intracellular proteins; and those with functions in pH 
maintenance in the retina, in metabolism, and as chaperones. 
The protein with by far the largest number of mutations is 
rhodopsin (>100 mutants), while the others range contain 
from 40 mutations (for the retina-specific crumbs homolog 1 
[CRB1]) to 1 (for transducin oc) (see Supplementary Material 3 
and Supplementary Figure S10). Structural analyses of the 
different mutations mapped on the available structures or 
homology models (156 mutations) indicated that the majority 
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of these are within the hydrophobic core of the corresponding 
proteins and are therefore likely to cause misfolding (see 
Supplementary Material 3). Mapping all proteins involved in 
vision-related diseases into the network made it apparent that 
the core visual pathway is the most susceptible to disease, and 
that the other functional modules are relatively robust. Out of 
36 proteins considered here to be involved in retinal 
degeneration, the majority (20 proteins) are localized in ROS 
(the remaining are found in other regions of the rod cells or in 
other cells involved in retina homeostasis (pH control), retinal 
recycling, or phagocytosis of the ROS discs). We found two 
cases of an ROS protein also expressed in other tissues, with no 
other apparent phenotype; e.g., isocitrate dehydrogenase 
NAD-dependent subunit B is found in many cell types besides 
rod cells. 

The prevalence of proteins from the core visual pathway in 
disease may have several explanations: first, mutations in other 
modules central for cellular function may result in a systemic 
all-or-nothing behavior, affecting the overall viability or proper 
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development of an organism and thereby causing early death. 
This may be true for critical cytoskeletal proteins and GTPases, 
as e.g., those involved in vesicle trafficking and maturation, and 
for proteins involved in metabolic activity. Second, the lack of 
redundancy for the very specific functions within the visual 
pathway might cause it to be more susceptible. Here, evolution 
may have favored high-end functional properties over the 
robustness of the pathway. Thus, lack of redundancy may have 
been accepted by evolution even though it interferes with 
robustness as a pay-off for the high-end performance that is 
achieved in photoreceptors with single photon detection and 
with multi-color vision. 



Conclusions 

Taken together, this work suggests that rhodopsin is able to 
trigger several distinct physiological activities in addition to its 
primary function of closing and opening the cAMP-gated 
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cation channel. Considering protein interactions as a result of 
domain interactions has allowed us to increase the resolution, 
define discrete functional modules, and add a spatial dimen- 
sion to this network. Based on this study, we obtained a novel 
biological insight that offers new testable hypotheses, which 
have been partially validated through the experiments 
performed here, namely, of the connectivity of rhodopsin to 
small GTPases involved in cytoskeleton assembly/disassem- 
bly and dynamics, and to vesicle and Golgi trafficking. This 
suggests a role for rhodopsin in self-regulating and fine-tuning 
the structural and functional integrity of photoreceptors. 
Cytoskeleton changes, such as microtubule assembly reorga- 
nization, are likely to affect protein transport between the 
inner and outer segments during light-to-dark changes 
(reviewed in Reidel et al, 2008), as well as to regulate cell 
polarity and disc development. The involvement of rhodopsin 
in regulating intracellular Ca 2 + levels suggests its role in an 
overarching Ca 2 + -dependent regulatory network that deter- 
mines dynamic changes in kinase activity and protein complex 
assembly. This in turn results in higher order physiological 
behavior, such as cytoskeletal dynamics and vesicular 
trafficking tuned by light. At a systems level, these network 
relationships imply a concerted regulation of outer segment 
structure, polarity, and vesicular trafficking orchestrated via 
GTPase-guided signaling pathways activated by light, Ca 2 + - 
regulated processes activated by cGMP-gated channel activity 
(and thus also by light), and cytoskeletal and ciliary dynamics 
(which may also be fine-tuned by light). With respect to 
disease, we can conclude that, among at least four pathways 
driven or regulated by rhodopsin, the visual pathway is the 
only one highly associated with disease, whereas all others are 
relatively unaffected. Conceptually, our work presents a 
general approach applicable to the analysis of any cellular 
pathway. The resulting comprehensive multiscale 'vision 
network' can serve as a basis for elucidating physiological 
principles of photoreceptor function and may help to identify 
potential disease-associated proteins and to guide signaling 
branch-specific therapy development. 

Materials and methods 
Isolation of ROS and ROS discs 

Porcine eyes were obtained from a local slaughterhouse. After the 
retinae were dissected, two approaches for ROS isolation were 
compared: that to Molday et al (1987) with that of Papermaster and 
Dreyer (1974). Briefly, for the Molday protocol, ROS was detached 
from the retinal tissue by gentle mechanical homogenization in cold 
isolation medium (20% (w/v) sucrose, 20 mM Tris, 2mM MgCl 2 , 
130 mM NaCl, at pH 7.2) and separated from the homogenate by 
loading onto a 27-50% linear sucrose density gradient. Alternatively, 
fresh retinae were homogenized by shaking in cold isolation medium 
(34% (w/v) sucrose, 65 mM NaCl, 2 mM MgCl 2 , and 5 mM Tris-acetate 
buffer, pH 7.4). ROS was then pelleted by centrifugation, and the 
remaining retinal tissue was rehomogenized with a teflon homo- 
genizer. Supernatants from both homogenization steps (crude ROS) 
were combined and loaded onto step-density gradients of 1.15, 1.13, 
and 1.11 g/ml sucrose. After cold centrifugation in a Beckman SW40 
rotor for 1 h at 38 000r.p.m., purified ROS was collected from the 
surface of a 1.11-1.13 g/ml sucrose gradient, and the protein content 
was determined by Bradford assay (Bio-Rad). Osmotically intact discs 
were isolated from ROS according to Smith et al (1975). ROS was 
ruptured by osmotic shock and intact discs were separated by flotation 
in 10% Ficoll (Sigma). After centrifugation (120 000 g, 2 h, 4°C), intact 



discs were harvested from the Ficoll surface. The purity of the ROS 
preparations was checked either optically by microscope (Figure 1A, 
inset) or by immunoblot analysis for RIS markers (BIP and Tom20) 
(Supplementary Figure Sll). 



Sucrose density gradient centrifugation 

ROS or intact discs were ruptured by osmotic shock, and the 
membranes were separated from soluble fraction by centrifugation. 
An amount of membrane equivalent to 1 mg protein was solubilized in 
1% (w/v) (3-dodecylmaltoside (DM) (Sigma) as described (Mueller 
and Eichacker, 1999), loaded onto linear 0.1-1.0 M sucrose gradients, 
and centrifuged for 17 h at 230 000 g at 4°C. Individual gradient 
fractions were either loaded directly for SDS-PAGE or were precipi- 
tated with methanol/chloroform as previously described (Wessel and 
Fliigge, 1984). 



SDS-PAGE and immunoblotting 

SDS-PAGE and subsequent immunoblotting onto PVDF membranes 
(Amersham) were carried according to standard procedures. Anti- 
body-antigen complexes were visualized using enhanced chemilumi- 
nescence detection (ECL + , Amersham) on Hyperfilm (Amersham) . 
Immunoblots were incubated with the following antibodies: anti-RhoA 
26C4 and anti-ROCK II H-85 (Santa Cruz), anti-visual arrestin, anti- 
transducin a, and anti-rhodopsin (Affinity BioReagents), anti-rhodop- 
sin (Acris Antibodies), anti-CRMP-2 (C4G, a generous gift from M 
Morishima and Y Ihara, University of Tokyo, Japan), anti-Racl (BD 
Transduction Laboratories), anti-RhoABC (Sigma), anti-BIP (BD 
Bioscience), and anti-Tom20 (BD Bioscience) . HRP-coupled secondary 
goat-anti-rabbit and goat-anti-mouse antibodies were obtained from 
Jackson ImmunoResearch. 



Immunoprecipitation 

IP was performed with anti-RhoA-agarose- or anti-Racl -agarose- 
conjugated antibodies (Santa Cruz) or anti-rhodopsin (Acris Anti- 
bodies). An amount of ROS equivalent to 500 |ig protein was ruptured 
by osmotic shock in lysis buffer (50 mM NaCl, 1 mM EDTA, 20 mM 
Tris-HCl, pH 6.8) and centrifuged to separate the membrane and 
soluble fractions. The membrane fraction was solubilized in 1 % (w/v) 
DM, and the soluble fraction was directly subjected to IP. For anti- 
rhodopsin IP, solubilized ROS (1% DM) was directly subjected to IP. 
Non-specific protein binding to agarose beads was prevented by pre- 
incubation of the fraction with 25% protein G-agarose (Santa Cruz). 
For IP, the ROS fractions were incubated with 5-10 (ig of antibody 
conjugate/antibody at 4°C for 3 h or overnight with rotation. As a 
control for non-specific antibody binding species, specific IgGs (Sigma- 
Aldrich) were used. 



Immunohistochemistry 

Porcine eyes were obtained from a local slaughterhouse, fixed in 4% 
paraformaldehyde in 0.1 M phosphate buffer (PB) for 4 h, and rinsed in 
0.1 M phosphate-buffered saline (PBS). Cornea, lens, and vitreous 
body were removed, and the retina was cut in 1 .5 x 1 .5 cm pieces. The 
fixed tissue was cryoprotected at 4°C stepwise in 10, 20, and 30% 
sucrose in PBS, for 1 h for the first two steps and overnight for the last 
step. Retina was then embedded in tissue-freezing medium (Leica 
Microsystems) and frozen in liquid nitrogen. In all, 12|im sections 
were prepared, mounted on Superfrost glass slides, and air dried at 
37°C. Retinal sections were rinsed in PBS and then non-specific 
binding sites were blocked with PBS containing 10% normal goat 
serum (NGS), 1% bovine serum albumin (BSA), and 0.3% Triton 
X-100 for 1 h at room temperature. Sections were incubated overnight 
at 4°C with the following primary antibodies diluted in blocking 
solution: rabbit anti-CRMP-2 (1:300; Abeam), rabbit anti-Racl (1:100; 
Sigma), rabbit anti-ROCK II (1:300; Abeam), rabbit anti-PDE65 (1:200; 
ABR), or mouse anti-rhodopsin (1:200; Millipore). Sections were then 
washed in PBS and incubated with the appropriate fluorescent-labeled 
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secondary antibody (goat-anti-rabbit IgG-Alexa 568 or goat anti-mouse 
IgG-Alexa 568; Molecular Probes) diluted 1:500 in PBS for 1 h at room 
temperature. Nuclei were counterstained with Sytox Green Nucleic 
Acid Stain (Molecular Probes). After three final washes in PBS, 
sections were mounted with Mowiol 4-88 (Polysciences) . As negative 
controls, the primary antibodies were also omitted; in these cases, no 
staining was observed. Stained cryostat sections were analyzed and 
scanned with a confocal laser scanning microscope (Zeiss LSM510 
META, Jena, Germany), using an argon laser at 488 nm and a He/Ne 
laser at 543 nm excitation with appropriate filter sets. Images were 
taken sequentially to assure that only one channel was detected at a 
time. The Sytox Green nuclear stain was allocated to the blue color 
channel for convenient viewing. Transmitted light images with DIC 
optics (Nomarski) were recorded simultaneously. Control sections 
without primary antibody incubation were scanned with the same 
laser and detection settings. 



Concanavalin A pull down 

For concanavalin A pull-down experiments, an amount of ROS 
equivalent to 300 (ig protein was ruptured by osmotic shock in 
lysis buffer. The membrane ROS fraction was solubilized in 1 % (w/v) 
P-DM (Sigma- Aldrich) , and the soluble fraction was directly subjected 
to concanavalin A pull down. The ROS fractions were incubated with 
50 ul of concanavalin A sepharose (Amersham Biosciences) conjugate 
for 3 h at 4°C. Non-specific protein binding of the rhodopsin-associated 
protein to concanavalin A was prevented by performing the pull down 
in the presence of 0.2 mM a-methylmannoside (as the presence of 
a-methylmannoside lowers the affinity of proteins for the beads) . 



Blue-native PAGE 

Membranes from either ROS or intact discs corresponding to 300 ug 
protein were suspended in 60 ul buffer containing 750 mM 
8-aminocapronic acid, 50 mM bis-Tris, pH 7.0, and 0.5 mM EDTA, 
and then solubilized in 1% (w/v) DM. The solubilized membrane 
samples were added to a buffer containing 5 % (w/v) Serva Blue G in 
750 mM 8-aminocapronic acid, loaded onto 4-12% PAA gradient gels, 
and electrophoresed (Schagger and von Jagow, 1991). To separate in 
the second dimension, gel lanes were incubated for 20min in 
solubilization buffer containing 2% (w/v) SDS, 66 mM DTT, and 
66 mM Na 2 C0 3 , and loaded onto denaturing PAA gels. 



MS-MALDI-TOF 

Selected spots were excised from dried silver-stained gels, destained 
(GharahdaghietaZ, 1999), dehydrated in 40 % acetonitrile (100 ul), and 
subjected to tryptic proteolysis in 1 mM Tris-HCl, pH 7.5, and 0.01 ug/ul 
trypsin. In parallel studies, proteins excised from dried gels were 
subjected to SDS removal by ion-pair extraction before in-gel tryptic 
proteolysis as described (Zischka et al, 2004) . MALDI-TOF PMFs were 
obtained on a Bruker Reflex III mass spectrometer (Bruker Daltonics, 
Bremen) . Aliquots from each tryptic digest were co-crystallized with a 
matrix composed of 2.5-dihydroxybenzoic acid (20mg/ml in 20% 
acetonitrile, 0.1% trifluoro acetic acid (TEA)) and 2-hydroxy-5-methox- 
ybenzoic acid (20mg/ml in 20% acetonitrile, 0.1% TFA) in a 9:1 ratio 
(v/v) on 400 um AnchorChip targets (Bruker Daltonics) . Alternatively, 
PMF and MS/MS spectra were measured on AB4700 mass spectrometer 
(Applied Biosystems, Darmstadt, Germany), and aliquots from each 
tryptic digest were co-crystallized with a matrix comprising 5% 
cyanohydroxycinnamic acid (in 70% acetonitrile, 0.1% TFA) on steel 
targets (Applied Biosystems). Database searches were performed using 
the Mascot software (Perkins etal, 1999). 



MS-Orbitrap 

LC-MS/MS analysis was performed on an Ultimate3000 nano-HPLC 
system (Dionex) coupled to an LTQ OrbitrapXL mass spectrometer 
(Thermo Fisher Scientific) by a nanospray ion source. Tryptic peptide 
mixtures were automatically injected and loaded with a flow rate of 



30 ul/min in 95% buffer C (0.5% TFA in HPLC-grade water) and 5% 
buffer B (98% acetonitrile, 0.1% formic acid in HPLC-grade water) 
onto a nanotrap column (100 um i.d. x 2 cm, packed with Acclaim 
PepMaplOO C18, 5 um, 100 A, Dionex). After 5min, peptides were 
eluted and separated on an analytical column (75 um i.d. x 15 cm, 
Acclaim PepMaplOO C18, 3 um, 100 A, Dionex) by a linear gradient 
from 5 to 40% of buffer B in buffer A (2% acetonitrile, 0.1% formic 
acid in HPLC-grade water) at a flow rate of 300 nl/min over 90 min. The 
remaining peptides were eluted by a short gradient of 40-100% buffer 
B over 5 min. Eluting peptides were analyzed by the LTQ OrbitrapXL. 
From the high-resolution MS pre-scan with a mass range of 300-1500, 
the 10 most intense peptide ions were selected for fragment analysis in 
the linear ion trap if they exceeded an intensity of at least 200 counts 
and if they were at least doubly charged. The normalized collision 
energy for CID was set to a value of 35, and the resulting fragments 
were detected with normal resolution in the linear ion trap. The lock 
mass option was activated, and a background signal of a mass of 
445.12002 was used for the lock mass. Every ion selected for 
fragmentation was excluded for 30 s by dynamic exclusion. The raw 
data were analyzed using Sequest (Thermo Fisher Scientific) and 
Scaffold (Proteome Software) as described previously (Gloeckner et al, 
2009) against a non-redundant pig, human, mouse, rat, and bovine 
protein sequence database derived in-house from UnireflOO, due to an 
insufficient number of entries for porcine proteins in the databases. 
Proteins were considered to be specific when they displayed two or 
more peptides (with a peptide probability > 95 %) in at least two out of 
four experiments. The protein probability threshold was set to 99%. 
Contaminants such as keratins were removed. 



Comparison of different proteomic data sets 
determined in ROS 

All proteins identified in the three different proteomic data sets were 
mapped to their corresponding human ortholog gene IDs by sequence 
comparison (using the default BLAST value of 10) and then compared. 
Since the proteomic analysis of Liu et al (2007) also contains part of the 
cilium, our 'near-to-complete' proteomic data set was defined as the 
union of the protein set identified by Kwok et al (2008) with the one 
determined here (Figure IB). 



Protein interaction network analysis 

All the results described in our studies were uploaded into 
Supplementary Tables S3 and S6 according to standard database 
curation rules (Zanzoni et al, 2002; Ceol et al, 2010) . Results from pull- 
down and co-IP experiments were resolved as binary protein 
interactions, in which each bait protein was linked to all identified 
preys. Co-sedimentation and complex-purification experiments that 
unambiguously identified protein complexes but lacked sufficient 
detail to determine their exact interaction topology are represented in 
the database as a list of interactors (complex members). A 
comprehensive literature mining and database curation effort was 
also carried out to include as exhaustively as possible the set of 
rhodopsin/vision-related interactions already described in the scien- 
tific literature. The curated interaction sets were represented and 
analyzed by the Cytoscape visualization and analysis software 
(Shannon et al, 2003). PPI data from databases included interactions 
determined from ROS extracts (by co-sedimentation or affinity 
chromatography). However, in the majority of cases, interaction 
information was derived from in-vitro experiments, such as large-scale 
yeast two-hybrid screens, or tandem-affinity purifications in artificial 
cell systems. Additionally, low-scale PPI data from the literature 
include data determined with quantitative affinity methods, such as 
isothermal titration calorimetry, surface plasmon resonance, nuclear 
magnetic resonance, and peptide arrays (using purified proteins), or 
PPI data from non-quantitative methods, such as affinity chromato- 
graphy (GST pull down), crosslinking, and enzyme assays. According 
to the MINT curation rules, interactions were considered to be direct if 
they were supported with evidence obtained with one of the following 
methods, as described in the PSI MI controlled vocabulary: two- 
hybrid, enzymatic studies, two-hybrid pooling approach, two-hybrid 
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array, p-lactamase complementation, surface plasmon resonance, 
fluorescence resonance energy transfer, biochemical, biophysical, 
protein arrays, protease assays, bimolecular fluorescence complemen- 
tation, far-western blotting, crosslinking studies, electron paramag- 
netic resonance, two-hybrid fragment pooling approach, protein 
kinase assay, GTPase assay, enzyme-linked immunosorbent assays, 
peptide arrays, isothermal titration calorimetry, bioluminescence 
resonance energy transfer, competition binding, fluorescence technol- 
ogies, antibody arrays, saturation binding, fluorescence polarization 
spectroscopy, protease accessibility laddering, affinity technologies, 
protein crosslinking with a bifunctional reagent, ubiquitin reconstruc- 
tion, fluorescence microscopy, p-galactosidase complementation, 
biochemical activity, classical fluorescence spectroscopy, fluorescence 
technology, phosphatase assay, and reconstituted complex. 



Structural information and interaction modeling 

Structural information was derived by a combined approach of 
comparing different domain interaction types, as listed in the 3 Did 
database (http://3did.irbbarcelona.org/; Stein et al 2005, 2009), 
between two interacting proteins for which there was experimental 
evidence that they could form a complex. The 3DID database was 
improved by analyzing all structures for crystallographic artifacts, 
using: (i) the interaction annotation of the author, or if this was not 
available and (ii) the protein quarternary structure (PQS) method 
(Henrick and Thornton, 1998). The confidence of two domains to 
mediate the interaction was then assessed using the InterPreTS 
(http://www.russelllab.org/cgi-bin/tools/interprets.pl/; Aloy and 
Russell, 2003) scoring system, which evaluated sequence similarity 
and amino-acid propensities in the interface. We further screened for 
all X-ray and NMR complex structures and homologs (with a sequence 
similarity threshold of 70 % ) among all 360 proteins of the network. For 
further details, see Supplementary Material 1. 



PDE5 subunit activity assay 

Recombinant PDE5 protein (rhPDE5) was obtained from GenWay 
Biotech at a concentration of 0.7 u.g/u.1 in storage buffer (10 mM Tris, 
pH 8.0, 0.1% Triton X-100, 0.002% NaN 3 , and lOmM dithiothreitol) . 
An amount of ROS (in isolation medium) corresponding to 100 ug 
protein was ruptured by three freeze-thaw cycles in liquid nitrogen and 
centrifuged at 4°C for 30min at 100 000 g (Beckman Optima 
ultracentrifuge; Rotor TLA110). The resulting pellet, containing the 
membranous fraction, was resuspended in 100 ul incubation buffer 
(25 mM Hepes, 20 mM Tris-HCl, pH 7.5, 1 mM dithiothreitol, 1 mM 
MgCl 2 , 5mM EDTA, 150 mM NaCl, and protease inhibitor cocktail 
(Roche)) and then incubated with different amounts (0, 0.5, 1,2,4, 6, 
or 8 (ig) of rhPDE5 for 1 h at 37°C in a horizontal shaker. To rule out that 
the Triton X-100 in the storage media affected the recombinant PDE5, 
all samples were adjusted to the same volume (volume of the sample 
with the highest PDE5 concentration used) with storage buffer. 
Samples were separated into membrane and soluble fractions by 
centrifugation at 4°C for 30min at 100 000 g and analyzed by SDS- 
PAGE and western blot using anti-Racl antibodies and anti-PDE5 
antibodies. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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